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The dynamical phase diagram of the fractional Langevin equation is investigated for harmonically 
bound particle. It is shown that critical exponents mark dynamical transitions in the behavior of 
the system. Four different critical exponents are found, (i) Qc = 0.402 ±0.002 marks a transition to 
a non-monotonic under-damped phase, (ii) an = 0.441... marks a transition to a resonance phase 
when an external oscillating field drives the system, (iii) a^i = 0.527... and (iv) a^^ = 0.707... marks 
transition to a double peak phase of the "loss" when such an oscillating field present. As a physical 
explanation we present a cage effect, where the medium induces an elastic type of friction. Phase 
diagrams describing over-damped, under-damped regimes, motion and resonances, show behaviors 
different from normal. 
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I. INTRODUCTION 

In this paper we investigate the phenomenological 
description of stochastic processes using a Fractional 
Langevin Equation (FLE) [TJ H E HH [7] . While in 
simple systems the memory friction kernel is an exponen- 
tially decaying function or a delta function, in complex 
out of equilibrium systems the picture is in some cases 
different. Namely the relaxation is of a power law type, 
and the particle may exhibit anomalous diffusion and re- 
laxation (Hj. Mathematically such systems are modeled 

using fractional calculus, e.g. An example is the 

recent experiment on protein dynamics of the group of 
Xie |9j [1^ . There anomalous dynamics of the coordinate 
X, describing donor-acceptor distance was recorded, and 
a FLE (see Sec. [n]) was found to describe the experimen- 
tal data. The motion of x is bounded by a harmonic force 
field, and the equation of motion for the average {x) is 



{x)+uj'^{x) -1-7 



(P{x) 
dt°' 







(1) 



where < a < 1, a; is the harmonic frequency and 7 > 0. 
For the case a = 1 we get the usual damped oscilla- 
tor [TT]. For such a normal case two types of behav- 
iors: the under-damped and the over-damped motions 
are found. In the under-damped case {x) is oscillating, 
and crossing the zero line, while for over-damped case 
[x] is monotonically decaying with no zero crossing. For 
a = 1, there exist a critical frequency = 2 ^^lich 
separates these two types of motion. Here we explore 
a similar scenario for the fractional oscillator, and find 
rich types of physical behaviors. It is known [T2] that 
in the long time limit all solutions {x) (i.e. any < a, 
< 7 and < w ) decay monotonically, some what like 
the over-damped behavior of the usual oscillator, how- 
ever now the decay is of power law type. The interesting 
Physics occurs at shorter times where the solution may 
exhibit different types of relaxation and oscillations. We 
find that for a < Uc the solution is non-monotonic for 
any set of parameters (7, w > 0). Thus we find a critical 
a which marks a dynamical transition in the behavior of 
the system. 



We also investigate the response of a system described 
by Eq. ([T]) to an external oscillating force Fq cos{U,t). For 
the regular case of a = 1, a resonance is present if the 
frequency uj is larger than the critical value ^jspl. The 
behavior for < a < 1 is quite different and we dis- 
cover that the transition between a — > 1 and a ^ is 
not smooth. In particular we find another critical expo- 
nent a^j, where for a < a resonance is always present. 
Other critical exponents are found for the imaginary part 
of complex susceptibility. Our goal is to clarify the na- 
ture of solution to Eq. (jlj, investigate the meaning of 
fractional critical frequency with and without external 
oscillating force and provide a mathematical tool box for 
finding and plotting solutions of Eq. ([T|). Our finding 
that critical as mark dynamical transitions is very sur- 
prising and couldn't be obtained without our mathemat- 
ical treatment. 

The paper is organized as follows. In Sec. ITT] we present 



the FLE. In Sec. Ill we present two different methods for 



solution for first order moments and examples are solved 
in Sec. |IV[ In Sec. |V] we present different definitions 
for over-damped and under-damped motion and find the 
critical exponent a,,- We also interpret our results from a 
more Physical point of view and discuss the cage effect as 
a viscoelastic property of the medium. In Sec |VI| we in- 
troduce external oscillating force into the system and find 
the response for such force for the free (Sec. VI A I and 
harmonically bounded (Sec. VI B) particles. In Sec. VIC 
we investigate the properties of the "loss" - the imaginary 
part of complex susceptibility. A summary is provided in 
Sec. |VII| and the three Appendixes deal with some tech- 
nical aspects. A brief summary of some of our results 
was published 



II. THE MODEL 



We consider the FLE 



d^x(t) 
m — V 7 



1 dx V 

^dt^F(x,t) + i{t) (2) 



{t - ty dt 
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where 7 > is a generalized friction constant (7 = 
^7r(l — a)), < a < 1 is the fractional exponent and 
£,{t) is a stationary, Fractional Gaussian noise [TU [TS] 
satisfying the fluctuation-dissipation relation |16j 



m) = o, mat)) = hT^\t - 1\~ 



(3) 



F{x,t) is an external force. We follow expe rim ent (TU] 



-mup'x, later in Sec. VI 



we will 



and assume F(x,t) 
have F{x, t) = —muj^x + Fq cos{ilt). In Laplace Space it 
is easy to show, using the convolution theorem 



x{s) 



1 



as) 



s2 + is/3(s)+w■ 



^«o 



s2+i-s/3(s)+C^2 

where xq and vq are initial conditions and 
f3{s) = ^T{1 - a)s°'-\ 



(4) 



(5) 



All along this work the variable in the parenthesis de- 
fines the space we are working in (e.g. x{s) is the 
Laplace Transform of x{t)). Eq. ^ with a = ^ and 
F{x, t) — —muPx, describes single protein dynamics [10]. 
An experimentally measured quantity is the normalized 
correlation function 



C^{t) 



{<t)xm 

(0.(0)2) ■ 



(6) 



In what follows, thermal initial conditions are assumed 
{i{t)x{Qi)) = 0, (x(0)2) = kbT/muj^ and (x(O)w(O)) = 0. 
From Eq. ^ we find 



4(s) 



- 7s 



a-l 



rygl- 



(7) 



It is easy to show that Cx{t) satisfies the following 
fractional Eq. 



dt°' 



= 0, 



(8) 



with the initial conditions Cx{0) = 1 and Cx{0) = 0, 
where the fractional derivative is defined in the Caputo 
sense [TTIIIH] 



(9) 



and o^T ^ is the Riemann-Liouville fractional opera- 
tor \mm 



oor'fit) 



r(i 



(i-t)-"/(i) dl 



Note that another way to write Eq. ([2| is 

m 



d^x , 

X + 7— h OJ^X ■ 

dt"' 



(10) 



(11) 



hence the name Fractional Langevin equation is justified. 
For the force free particle {F{x,t) = in Eq. Q) {x"^) cx 
t" [BJIin], this is a sub-diffusive behavior since < a < 1. 

The FLE in general and Eq. ([s]) in particular can be de- 
rived from the Kac-Zwanzig model of a Brownian particle 
coupled to a Harmonic bath |20] . Eq. (|8| can be derived 
also from the Fractional Kramers Equation {TP . The 
use of fractional-differential equations like Eqs. ( 8|11 1 be- 
came quite common in recent years [5J especially 
in the context of anomalous diffusion [8j ESj El] . Sev- 
eral other fractional oscillator equations were considered 
in the literature [251 ESI EZ] and general solutions for 
fractional-differential equations of the type Eq. J8 ) are 
studied in the mathematical literature [HI EH [29 . In 
the next section we present a practical recipe, a tool box, 
for the solution of Eq. ([s]), and show how to plot its 
solution. Our methods are general and can be applied 
to other linear fractional differential equations. From a 
more physical point of view the following questions are 
addressed in the next sections: (i) When is Cx{t) pos- 
itive? (i.e. similar to over-damped behavior, a = 1). 
(ii) When is Cx(t) non-monotonic? (similar to under- 
damped motion for a = 1) (iii) What is the analogue to 
the critical frequency lOc = ^ found for the a = 1 case, 
(iv) Does a critical exponent ac exist and if so, what is 
it's value? 



III. THE GENERAL SOLUTION 

In this section we will present a recipe for an analytical 
solution of Eq. (|8|. 

The formula for Cx{t) in Laplace space is Eq. ([t]) 
and if a was an integer then performing an Inverse 
Laplace Transform would be simple using analysis of 
poles [SnilSI], because then the denominator of Eq. Q is 
a polynomial. We assume a is of the form a = -, where 
q > p > are integers and ^ is irreducible (i.e not equal 
to some other - where I < p and n < q are integers) . 



Method A 



We rewrite Eq. ([t]) as 

Q{s) 



Cx{s) 



- -fS 1 



s2 + 759 + LU^ Qis) 



(s + 7s? ^) Q(s) 



Pis) 



(12) 

where P{s) is a polynomial in s. According to a mathe- 
matical theorem |18j we can always find Q{s) such that 
the denominator of Eq. ( 12 ) 



P(s) = (s^ +-fs-^ +uj^)Q{s) 



(13) 



is a regular polynomial in s. The polynomial Q{s) is 
called the complementary polynomial. The task of find- 
ing Q{s) is simple, Q{s) is a polynomial in of degree 



3 



2qiq - 1) 



Q{s) 



29(9-1) 



we can write Eg. (plj ) by the means of generalized Mittag- 
LefHer function 1331 



(14) 



with 529(9-1) = 1- The coefficients are found by 
equating all the coefficients in the expansion of the prod- 
uct (5(s)(s^ + 759 with non-integer powers of s, to 
zero. This produces a linear system of 2q(q—l) equations 
for 2q(q ~ 1) variables, which in principle is a solvable 
problem. 

We also assume that all 2q zeros of P{s) are distinct. 
The generalization to the case where the zeros are not 
distinct will be treated later, such a behavior is related 
to a critical frequency of the system. We use the partial 
fraction expansion 



1 

-r, *^^i,n 

The Mittag-LefELer function satisfies 



y" 



with 



r(M - v) 



00. 



(22) 



(23) 



(24) 



P{s) 



2q 

E 



A, 



(15) 



where a}, are the solutions of P{s) = and are con- 
stants defined as 



Actually we have finished our task, we can now perform 
an Inverse Laplace Transform of expressions like Eq. ([t]) 
and even more general expressions which can be pre- 
sented as a fraction of polynomials with fractional pow- 
ers. To summarize, we need to follow four steps (i) Cal- 
culate Q{s), which is equivalent to the diagonalisation 
of a matrix of size {2q{q — 1)) . (ii) Find the zeros of 
P{s), Ofc, and th e co efficients of partial fractions expan- 



Ak = 



1 



dP{s) 



(16) 



ds 



sion ,Ak of Eq. (|16|. (iii) Find the coefficients Bmj for 



It can be shown [15" that 
29 



Eq. ( 18 1 and write Cx{s) as the sum in Eq. ( 19 1, (iv) Use 
Eqs. ( 20|21|22" ) to inverse Laplace transform Eq. ( |l9| ), we 
find 



P{s) 



E 

k=l 



/I™ A , 



m ~ 0, 1, 



ak 



,2q~l. 



(17) 



29-1 9-1 29 m-j/q 

^^w = EEE^ 

m=0 j=0 k=l 



R ■ A, 



r(^) 



M-.dkt) (25) 

q 



The numerator of Eq. ( 12 1 is written using the expansion where for j = 

29-1 9-1 

g(,) (s + 7st-i) = yj^,s"^-i. (18) 



7(-, akt) 

' , = 1, or using Eq. (|22| 

\q) 



m=0 j=0 



Hence using Eqs. ( |12|15|18 l 

29-1 9-1 29 



29-1 9-1 29 

m=0 j=0 k = l 



BmjAkti E, , , j (afct) 

' ' a 



(26) 



c,{s) = y Y^Y^<Bm,Ak 



S 1 



m=Q j=0 k=l 



ak 



(19) 



So finally we reduced the problem of calculating the In- 
verse Laplace Transform of Eq. ([T]) to performing Inverse 
Laplace Transform for 



and 



1 



si{s - ak) 



(20) 



— ^7(-,M) (21) 



where j{'^,akt) is a tabulated Incomplete Gamma Func- 
tion [^32]. Using the series expansion for "f{a,x) 



and for j ~ E, , , j (a/ji) ~ e"'-'*. Now we have a prac- 

tical tool for finding an explicit analytical solution for 
the fractional damped harmonic oscillator. Since 7(0,0;) 
is tabulated in programs like Mathematica, the solution 
which is a finite sum of such incomplete gamma func- 
tions, can be plotted. 



B. Method B 

The task of finding Q{s) is sometimes difficult since as 
described in Sec. |III A| generally one must solve a linear 
system of 2q(q — 1) equations with 2q(q — 1) variables. 
But for the special case of Eq. (j?]) we will provide a simple 
method for finding Q{s) and P{s). By writing 



a(t) 



S + JSl' 



(s2+Cj2)^ + (-l)9-l79sP 



(s2 + Cj2)9 + (-1)9-1795^ 



n=0 



T{a + n+l) 



+ 7s 9 + OJ^ 



(27) 
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1 



FIG. 1: The short time behavior for Cj:{t) with o = | and 
7 = 1, versus t. Three types of solutions are presented (a) 
u = 0.3 and the function decays monotonically. (b) u = 
LOz ~ 1.053 the transition between motion with and without 
zero crossing, Cx{t) = at a single point in time, and Cx{t) 
does not cross the zero line, (c) uj — 3 the under-damped 
regime. Note that for large t, Cxit) > 0. 



we can write 



Qis) 



and 



(28) 



(29) 



One easily sees that indeed P(s) is a polynomial in s. 
As for Q{s), it is found by standard method of division 
of two polynomials in s^^' and it is also a polynomial in 



FIG. 2: The long time behavior for Cx{t) with a = | and 

7 = l,versus t, on a loglog scale. Three types of solutions are 

presented (a) u) = 0.3 over-damped limit, monotonic decay. 

(b) uo = LOz ~ 1.053 the transition between motion with and 

without zero crossing, (c) u = ?, the under damped regime, 

for large t, Cx{t) > (for short time Cx{t) < 0, not shown). 

Notice that non-monotonic decay in the under-damped case 

is observed only for short times, while for long-times Cx{t) oc 
1 

t~2. 



The first step is to find the complementary polynomial 
of the denominator on the right hand-side of Eq. (31 1. 



Using method A of the previous section one can easily 
see that the coefficients of the complementary polynomial 
Q{s) are 



and 



Since the fraction of two polynomials 



Bi = 1, B3 = S2 = 0, B, = -7, Bo 



Q{s) _ .^53 _^^2 



(32) 



(2/29-fcj2^'' + (-l)9-l79y9 



y 



2q 



-jyP 



(30) 



is also a polynomial in y, we see that any solution of 
^2 = is also a solution of (y'^'^ + w'^V + 



^^lyip — 0, and by performing a substitute y = 



s? into Eq. (30 1 we find that indeed Q{s) is a polynomial 
in s^/^. After finding Q{s) and P{s) explicitly we can 
return to method A and use Eqs. p^p6 1 in order to 
write down the final solution. 



IV. EXAMPLES a 



AND Q 



The a = \ case was recently measured in experiment 
[To] , and so it will be our first illustration for the method. 
From Eq. ([?]) we get 



S + ^S 2 

§2 -1-753 -|- ii;2 



(31) 



smce 
P 



{s) = Q{s) ( 



+ -fs'i + u?) = -I- 21^3^ - 7^s -I- w*. 



We rewrite Eq. (31 1 as 



(33) 



(^^S) = , ..2\2 — (34) 



(§2 -I- tj2)2 _ ^2g 



and notice that the degree of the denominator of Eq. ( 34 1 



i.e. P{s), is 4. Its zeros are easily found using Ferrari 
formula [32], and we call them a^j, k = 1,...,4. The 
coefficients of the partial fraction expansion Ak are given 
by Eq. ([161 



A, 



1 

Aak{al+uj'^) -72 



(35) 



The partial fraction expansion is found using Eq. ( 19 1 

3 14 



^^(^)-EEE 



m=0 j"=0 k=l 



S - ak 



S " 



(36) 
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and the i?„ij in Eq. (18 1 are found using the numerator and other elements of the matrix Bmj are equal 0. Usin, 



of Eq. (341 



Eq. ( 25 ) the solution is 



B: 



30 



1, Bio^uj , Bqq = -7 , Boi^ju! 



(37) 



k = l 



(38) 



r 




FIG. 3: The short time behavior for Cx{t) with a — j and 
7 = 1, versus t. Three types of solutions are presented (a) 
LJ = 0.3 and the function decays monotonically. (b) u = 
LJz ~ 0.965 the transition between motion with and without 
zero crossing, Cx{t) does not cross the zero line, (c) uj — 3 
the under-damped regime. 



The asymptotic expansions Eq. ( 39 1 provides the be- 
havior for long (and short) times, but the intermediate 
behavior is not obvious. Using the exact solution Eq. ( [38| 
we plot Cx{t) for various values of 7 and u; in Fig. [ijand 
Fig. [2] Three types of behaviors exist (i) Monotonic 
decay of the solution - Fig. [T]^a).(ii) Non-monotonic de- 
cay in the non-negative half of the plane, Cx{t) > 0, - 
Fig. [ijb). (iii) Oscillations of the solution, where Cx{t) 
also takes negative values - Fig. [TJc). These are typ- 
ical behaviors of the solution which we found also in 
other parameter set (not shown). From Fig. [l] we identify 
uj — = 1.053 as a fractional critical point, in the sense 
that if w > Uz we have zero crossings for Cx(t). For 
a = I there exist also another fractional critical point 
ujm = 0.426 where for oj < ujm Cx{t) is monotonically 
decaying - Fig. [T|a). We will soon discuss lJz and ujjn 
more generally. 



By using a series expansion and some algebra we find 
Ofc and Eq. ( 38 ) has the following asymptotic behavior 



In this subsection we demonstrate the solution for a = 
using method B. Our goal is to invert Eq. nil. 



Cx{t) 



7 



'JUJ 

rXi)' 



t5 



t 



a;2r(i) 



2r(i) 







00 



(39) 

We see that Cx{t) for long times decays as a power-law, 
which is the signature of slow relaxation and anomalous 
diffusion. The same asymptotic results are found by ap- 
plying Tauberian Theorems [31] to Eq. (31). 



Cxis) 



s + 7s 4 



(40) 



We find the complementary polynomial Q{s) and P{s), 
using Eq. (l28| and Eq. (l29| 



and 



4„3 



(41) 



4„3 



7 s 



-f 7S4 + Cj2 



=s^ - 7s 4 + Slu^s'^ + 7 s 2 - 2cj^7s -i 7''s ■» + Sto'^s^ + lu^^s^ - juj'^s « -I- (42) 



and so we can write Eq. (40 1 as 



Cx{s) = — -4 — . 43 
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FIG. 4: The short time behavior for Cx{t) with a = | and 
7 = 1, versus t. The three types of solutions are presented 
(a) uj — 0.3 and the function oscillates above zero, (b) uj = 
LJz = 1.035 the transition between motion with and without 
zero crossing, Cx(t) does not cross the zero line, (c) uj = ?> 
oscillations with zero crossing for short times, for long times 
Cx{t) > and no oscillations are observed. For this case 
solution with monotonic decay are not found since q = p < 
= 0.402. 



The degree of the denominator of Eq. (43) is 8, so the 
zeros of the polynomial P{s) could only be found numeri- 
cally using a program like Mathematica. As in the previ- 
ous subsection we call the zeros of P{s) ak , k = 1, . . . ,8 
and the coefficients of partial fraction expansion are 



found using Eq. ( 16 1 

Ak = - 



(44) 



Writing the partial fraction expansion using Eq. ( 19 1 



m=Q j=Q k=l 



k B m j Ak 



(45) 



where the coefficients B„i j are found using the numerator 
of Eq. (pl 



B = 



/ 





7^6 








\ 

























7^u;^ 






3uj^ 


























































\ 


1 











/ 



Using Eq. ( 25 I the solution is 



(46) 



fc=i 

f (»=<.^' +2.«4 +"^<.f ) H- (^'4 +^''4) + ^^fcfM' 



or using Eq. ( 26 1 



k=\ 

+Ak (jtiE^^5{akt) {lu^ + 2Lu^al + Lo^al) + jh^ E^ s^akt) (uj^ak+io^al) + Lu^altij^E^ 7{akt)j 



(47) 



(48) 



The behavior described by Eq. (48 1 is plotted at Fig. |3] similar to the behavior for the a = \ case Fig. [T] The 
The three typical types of behavior are shown, which are 
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values of critical points for a = | are uJm ~ 0.707 (and so 
for uj = 0.3 we observe in Fig. p[a) a monotonic decay) 
and ujz ~ 0.965, a case plotted in Fig. |3jb). Finally the 
case a = I is shown in Fig. |4j The difference between 
q; = I and the former cases is that for Q; = | uj„i — and 
so we never observe a monotonic decay of the solution. 
More generally, this behavior is obtained for any a < 
ac — 0.402, as we shall soon show. 



DEFINITION OF OVER AND UNDER 
DAMPED MOTION 



7 



u;2r(l-a) 



(50) 



where the large t expression was obtained in [TOl [HI US] . 

The applicability of Eq. ( 50 1 is possible only under two 



conditions, the first one is obvious from Eq. ( |50| and is 

,2\ 1/" 



t > 1, 



(51) 



while the second is 



As mentioned in the introduction, when dealing with 
the normal damped oscillator one gets two types of so- 
lutions - over-damped and under-damped, the transition 
between these two behaviors happens at some point lOc 
called the critical point. For the over-damped motion 
(x) > for any t when {x{t = 0)) > 0, and there arc no 
oscillations, and for the under-damped case {x{t)) oscil- 
lates and crosses the zero line. For the fractional oscil- 
lator, we notice that there are different types of behav- 
iors. From Figs. ( 1|3|4| one notices that for short times 
there is an oscillating behavior either with and without 
zero-crossings or a monotonic decay type of behavior. 
So as in the regular damped oscillator we need to de- 
fine the transition between these behaviors. We propose 
three definitions for the point of transition between over- 
damped and under-damped motions, these are based on 
the various definitions that exist for the regular damped 
oscillator and give the same result for a — 1. The first 
option is to take the frequency coc for which two solu- 
tions of -P(s) = coincide, i.e. an appearance of a pole 
of a second order for Cx{s), and the general solution of 
Sec. Ill must be modified, as explained in Appendix | VIII 



\a\t > 1, 



(52) 



The second option is to take the minimal frequency lUz 
at which the solution Cx(t) crosses the zero line and the 

third is to take the minimal frequency w,„ at which '^'~^2t^^ 
crosses the zero line (i.e. Cx{t) is no longer a monoton- 
ically decaying function). For regular damped oscillator 
LJc — ujz = ^m, but in fractional case this is generally not 
true. 

Another difference between the fractional oscillator 
and the regular one is the distinction between short and 
long time behavior when < a < 1. The asymptotic be- 
havior for general a is obtained using general properties 
of polynomial solutions ^18j with Eq. ( 23 1 and Eq. ( 24 ) 
or using the Tauberian theorem 



Cx{t) 



2 



r(5 - a) 



t (49) 



where |al is an absolute value of a zero of the enumerator 
in Eq. ([7| (s^ + 7s" -|- w^), it also could be found from 
P{s) — 0. Mathematically the second condition is the 
magnitude of the radius of convergence for the power se- 
ries expansion of Cx{s) near s = (Tauberian Theorem), 
and is given by the non- fictitious poles of Cx{s). From 
Eq. ( 50 1 , Cx{t) > and for large t decays as a power-law 
(see Fig.|2]). 

From a more physical point of view, the FLE formalism 
was used [ID] to describe the fiuctuation of the distance 
between a fluorescein-tyrosine pair within a single protein 
on a time scales of 1 msec up to 10^ sec, and a power- 
law decay was observed (lately a theoretical model of 
Fractons was proposed in order to explain such phenom- 
ena |40j). Recent molecular dynamics simulations |41j 
studied the fluctuations of a donor-acceptor distance for 
a single protein, for short time scales 10~^ sec and os- 
cillations of the autocorrelation function were observed. 
The scenario of oscillations for short-times and a power- 
law decay for long-times set's well with the description 
by our solutions of the FLE. 



A. Critical point uoc 



For the case a = ^ we have four solutions for P{s) — 
{au k = 1,...,4) which are plotted in Fig. |5] At 
one point two solutions coincide, we call this point the 
critical point at which lo ^ lOc- For oj — lJc, the solution 
for Cx(t) Eq. (38 1 m ust b e modified because our general 



method derived Sec 



III 



is not valid. For a = ^ and 



1 
2 

^ - c^c = 53W\/(3/2) (1/2) 1/37'/' (see Eq. |59| Eq. 
is (sec Appendix, Eq. (95l) 



(381 



m=0 j=0 



■ 3 
.k=l 



1,H 



, {akt) + a^A ( (t + ma^ ^)i^i,i+i {a^t) 



(53) 
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where oi j- 
by Eqjssl 



Eq. (35 



in Eq. (88 I 



0,2 7^ 13 and a4 = 03 are zeros of Pis) given 



Brnj are found using Ec 



,yl3 = -(^1 + A2) by Eq. (87| and A is defined 



^1 and A2 by 



We emphasize that the critical point Wc does not al- 
ways exist. In Fig. [o] we plot the 10 solutions of -P(s) = 
for a = |, and demonstrate that no two solutions coin- 
cide. We will soon show that for any odd q and even p 
the critical point lOc does not exist. 

Mathematically at critical points one of the coefficients 
of partial fractions expansion, in Eq. (16 1, diverges. 



This happens because when two ak coincide, -P(s) can be 
written as (s — at) G{s), where G{s) is some polynomial 
in s and G{ak) 7^ 0, and 



A, 



2{s - ak)G{s) + {s~ ak) 



2 dG{s) 
ds 



(54) 



diverges. And so in order to find such critical point two 
conditions must be satisfied 



P{s) = 0, 



dP{s) 
ds 



= 0. 



Using Eq. (29 1 with oe = '-^ and lu ^ ujc, 



2qs {s 



p(-l)«-i7''sf-i =0. 



Solving these Eqs. we find 



± 



and LOc must satisfy 



1 



25^ 



(2 — a)Q!2-"75 



(55) 



(56) 



(57) 



(58) 



(59) 



Eq. ( 58 1 and Eq. ( 59 ) are only valid for even q or even 
{q + p) {recall ct = ^ a nd | is irreducible), since for any 
other case Eqs. ( 56|57 l has no solutions. The -I- sign in 
Eq. ( 58 1 is for the case of even q and odd p and the — sign 



for odd q and p. To see this, insert the solution Eq. (58) 
in Eq. (56) and then we must have (—1)^^^ < and 
hence q even, since ^ is irreducible p is odd. Similarly 
for the — solution in Eq. ( 58 1 . From this discussion it 



becomes clear why there are no critical frequency ujc for 
a=| (Fig.g. 



B. Critical points uiz and ui,n 

We divide the phase space to three different regions. 
{I) < uj < ujm the region where Gx{t) decays mono- 
tonically (//) LUm < uj < lOz the region of non-monotonic 




-D B -0 4 -0.2 



0.2 0.4 OB 



FIG. 5: The four solutions of P{ak) = 
plane for a — |, 7 = 1 and < ui < 5. 
ui — ujc ~ 0.6873 we have as = 04. 



- in the imaginary 
At the critical point 



« 




-O.B -0.4 -0.2 D.2 0.4 O.B O.E 
Real(ai,) 



FIG. 6: The ten solutions of P{ak) ~ for a = |, and 7 = 1, 
in the imaginary plane, for < it; < 5. All ten solutions are 
different and do not coincide, namely for this case there does 
not exist a critical frequency uoc since g = 5 is odd. The origin 
corresponds to u = Q. 



decay while Cx{t) always stays positive (///) < w 
the region of non-monotonic decay while part of the time 
Gx{t) is negative. Similar to Eq. (59) we find from di- 
mensional analysis 



and 



uiz = K^(a)72- 



(60) 



(61) 



where Kz{a) and n^^a) depend only on a. By inves- 
tigating analytical solution Eq. ( p6| for various a and 
7 = 1 we obtain functions k.z{ol) and Km{oL). The result- 
ing phase diagram is presented in Fig. [7] One can readily 
see that lOz, uJm and uJc defined by Eq. ( |59[ ) all coincide 
only for the normal case a = 1. 

An interesting behavior is observed for Km(a), as can 
be seen in Fig.|7]a sort of phase transition occurs around 
a ~ 0.4. We used the general method developed in 
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Non-Monotonic 
decay without 
Zem-Cmssing 



0.1 D.2 D,3 OA 0,5 0.6 0.7 0.8 D,9 



FIG. 7: The phase diagram of the fractional oscillator. Phase 

(a) monotonic decay of the correlation function Cx (t) , phase 

(b) non-monotonic decay without zero-crossing and (c) oscil- 
lations with zero crossings. The boundary between (b) and 

(c) is ujz = Kz(a) (solid line -I- squares), the boundary be- 
tween (a) and (b) is tOm = Km (a) (solid line + circles). For 
a < ac — 0.402, the phase of monotonic decay disappears, 
namely we do not find over damped behavior. The dotted 
curve is the critical line given by Eq. ( 59 1 . All the curves 

2 = O^m = 7/2. 



are calculated for 7 = 1. For a = 1, cuc 



Sec. Ill and explored the behavior of 
in Laplace space is given by 



1 dCx(t) 
w'^ dt ' 



sC^is) - 1 



1 



7s" 



which 



(62) 



For a < 0.4 we always observed zero crossings even if we 
decreased w to 10"'' (7 = 1), while for a > 0.404 there 
were no zero crossings beneath some finite w > 0, as is 
shown in Fig. [7] So we can conclude that there exist a 
critical a, ac ~ 0.402 ± 0.002, where for a < ac Cx{t) 
does not decay monotonically even if the frequency of 
the binding Harmonic field a; — 0. Note that the phase 
diagram Fig. [7] also exhibits some expected behaviors: 
as we increase u we find a critical line above which the 
solutions are non-monotonic and exhibit zero crossing, a 
line represented by Kz{a). 

To find accurate values of ac we also used a method 
based on Bernstein theorem [34]. According to the the- 
orem if and only if f{t) is positive then for any integer 



o<(-iy 



(s > 0) 



(63) 



for any < a < 1. But for n = 2 



(-1) 



ds" (s2 + S" + ^2)3 

(6s^ + {9a - - 2)s" + (a^ + a)s^°''^ 
-L0^{2 + a{a-l))s°'-^) , 



(64) 



and in the limit 



one can easily show that for 
a ^ 0.071 Eq.(64) has negative values. Actually for a — 



0.01 it is easily verified (by plotting Eq.(64)) that for 
any a; > Eq. (64) would have negative values. So for 

v(2) 



2 we have a upper bound for a, a^ 



0.071, where 



for any a < a 



(2) dCx(t) 



dt 



crosses the zero line and the 



relaxation is non-monotonic. If one wishes to increase 
the accuracy for such an upper bound one should inspect 

the behavior of (—1)" ^{i^'' for higher values of n, for any 
n we define such a bound as a'^"\ Using Mathematica 
we can proceed to higher values of n and find the upper 
bound a^"^ for different uj. In Fig.jsjwe plotted the upper 
bounds as a function of n for various w, we see that as n 
grows the upper bound a^"-* converges to some value < 1. 
The results achieved by this method converge to values 
very close to the values displayed in Fig. [7] for example 
for a; = and n — 150 a^^^°^ — 0.394, compared with 
ac ~ 0.402, obtained by inspection of exact solution. 

A physical explanation for this interesting result is 
based on the cage effect. For small a the friction force 
induced by the medium is not just slowing down the par- 
ticle but also causing the particle a rattling motion. To 
see this consider the FLE Eq. (Ill in the limit a ^ 



mx 



+ mj{x — xo) + muP'x w ^(i) 



(65) 



where 0:0 is the initial condition. 



Eq.(65l describes har- 
limit yields 



monic motion and the friction 7 in this a 
an elastic harmonic force. In this sense the medium is 
binding the particle preventing diffusion but forcing os- 
cillations. In the opposite limit of a — > 1 



mx + mjx + muj^x ~ £,{t), 



(66) 



the usual damped oscillator with noise is found. So from 
Eq. (651 an oscillating behavior is expected, even when 
Lu —^ which can be explained by the rattling motion of a 
particle in the cage formed by the surrounding particles. 
This behavior manifests in the non-monotonic oscillating 
solution we have found for small a. Our findings that ac 
marks a non smooth transition between normal friction 
a ^ 1 and elastic friction a — )■ is certainly a surprising 
result. 



where f{s) is the Laplace pair of f{t). As in the previous 
paragraph in order to check the monotonicity of Cx{t), 
we will inspect -\ '^'"X^*^ for zero crossings, or by speak- 



ing in the language of Bernstein theore ni b y exploring 
(^_iyn d JJji) ^ —f{s) is given by Eq. (62l. Using the 

, and so it is 



scaling relation of Eq. (61 1 we can set 7 

and n 



easily checked that for n 



i(-i)"^>o, 



VI. RESPONSE TO AN EXTERNAL FIELD 

In this section we will explore the response of the FLE 
to an external time dependent force. The response of a 
system to an oscillating time dependent field naturally 
leads to the phenomena of resonances, when the fre- 
quency of the external field matches a natural frequency 
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n (iterations) 



a^C,441 




' ft - 0,35 
■ft = 0,l 



1 1 , 1 1 ' I 

0.5 1 1.5 2 2 5 3 3 5 

Q 



FIG. 8: Using the Bernstein Method we find the upper bound 
a'"' as a function of n, the number of the derivatives. We 
consider three different uj (uj = 0, uj = | and w = j^). The 
dashed lines present the values found numerically in Fig. [7] 
As n grows a convergence of the bound is achieved as is seen 
for = as q'"' converges to 5^ 0.402 ± 0.002. 



FIG. 9: The response of FLE to an oscillating field for a 
free particle and different a'a. For a < ur a resonance is 
observed. All the curves are plotted with 7 = 1. 



Unbounded Particle 



of the system. The response of sub-diffusing systems to 
such time dependent field was the subject of intensive 
research [351 IMl ISZ] • In particular fractional approach to 
sub-diffusion naturally leads to anomalous response func- 
tions commonly found in many systems e.g. the Cole- 
Cole relaxation [Ml |3H1 [33] ■ In the next three subsections 
we will explore the response of FLE with and without the 
Harmonic potential and also investigate the behavior of 
the imaginary part of the complex susceptibility, i.e. the 
dielectric loss, for these cases. 

Our starting point is Eq. ^ with F{x,t) — 
Fq cos(f2t) — rauP'x, performing an average we obtain 



(a;) +7 



The solution in the long time regime is 
Fo 



cos (fit). 



(x(t)) 



cos{n{t - t))h{t) dt, 



(67) 



(68) 



where h{t) is soon defined. Eq. (68) could be written as 

{x{t)) ^ R{n)cos{nt + 0{n)) t^oo. (69) 

The response R{il) and the phase shift 9{il.) are obtained 
by the means of complex susceptibility 



(70) 



where x'(^) ^'^d x"(^) ^re the complex and the 
imaginary parts of the susceptibility, respectively, and 
h{-in) = e'^^h{t) dt [Uj. For the response 



R{n)^\xm\ 



and 

9{^) — arctan 
for the phase shift. 



For the unbounded particle we set w = in Eq. (67 1 
and using Eq. ^ we obtain for x(^) 



x{^) = h{-in) 



(73) 



By the use of Eq. ( [Tl] ) we can explicitly obtain now the 
behavior of i?(ri) for any < a < 1. For the normal dif- 
fusion case when a — 1, the response i?(f2) is a decaying 
function of and no resonance is observed, but the pic- 
ture is quite different for < a < 1. In this sub-diffusive 
part the response R(Cl) is not always a monotonically 
decaying function and could obtain a maximum, i.e. a 
resonance, even for such a free motion. As is seen in 
Fig.|9]for small enough a, R{fl) has a resonance, we will 
show that the existence of the resonance does not depend 
on any other parameter but a. 

Writing down the response R{il) explicitly 



R{n) 



^4 _^ ^2^2a _ 27^22+" COS 



(74) 



We are looking for the solutions of dR{il.)/dil. — 0, and 
hence 

An% + 2aj^n^j^-'^-2{a + 2)jn%+^cos(^'^) =0, (75) 



so the solution is 

7 



(^1) = 25 ((« + 2) COS (^) ± ^{a + 2)2 cos2 {^) - 8a) , 

(76) 

and is the frequency for which the resonance is found. 
(72) When the discriminant on the right hand of Eq. ( |76| is 
greater than zero there will be a resonance . The discrim- 
inant has no dependence on 7, and is always positive for 
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(a) (b) 




FIG. 10: The response of FLE to an oscillating field for a 
harmonically bounded particle and different as. (a) a = 0.2, 
ci; = 1 and 7 = 1. (b) a = 0.2, u) = \ and 7 = 10. (c) a = 0.7, 
ci; = 1 and 7 = 10. (d) a — 0.7, lu — 7 and 7 = 10. 



a < aR = 0.441021. 
tion 



which satisfies the following rela- 



{aR + 2ycos'(^)-8aR = (77) 

So no resonance is found for a > aji and for a < an 
there is always exist an ^Ir > for which the response 
will exhibit a resonance. This result of a resonance for a 
free particle FLE is highly unexpected but sets well with 
our description of the friction force for small a's as an 
elastic force due to the cage effect. 



FIG. 11: Phase diagram of the response of the system with 
a harmonic field to a oscillating time dependent force field. 
Two simple behaviors are found either a resonance exists 
("Resonance Phase") or not ( "No- Resonance Phase"). For 
a < an = 0.441... a resonance exists for any binding har- 
monic field and any friction 7. 



Eqs. d7T]78 l 

^{UJ^ - 172)2 + ^2fl2a + 2(cj2 _ n^)jn" COS (^) 

(79) 

and using dR{il)/dil — 0, we find 



2.(^)^+[2.(4-l)-4]cos(-)^ 

-4te-l)-0. 



B. Harmonically Bounded Particle 

Now we treat the response function of the FLE with 
a harmonic field, i.e. of a fractional oscillator. Starting 
with Eq. ( 67 1 we set the initial conditions xq — vq — 
and in the long time limit < — *■ 00 we obtain again (a:) = 
R{n) cos(m -I- 6I(»)) where R{n) and 9{n) are defined by 



Eqs. ( 71|72 1. The complex susceptibility for such case is 
obtained 



1 



(78) 



Eq. (78) was already obtained 



m m] for the frac- 
tional Klein-Kramers equation [5T] in the high damping 
limit and is called a generalized Rocard equation |45[ |46] . 
For a = 1 we find the complex susceptibility of a normal 
damped oscillator. 

We are interested in the resonance points for the re- 
sponse to the applied field, i.e. points of maximum of 
R{fl), which generally depend on fi, 7 and lu (see Fig. 10 1. 
For the normal oscillator there is a resonance if the condi- 
tion lli > ^7 is satisfied. If this condition is not satisfied, 

i?(ri) is a monotonically decreasing function of fl. From 



The exploration of Eq. (|80| (see Appendix IX I leads to 
two findings which are presented in Fig. |ll[ The first one 
is the existence of the same critical aR given by Eq. ( 77 1 
for the response of FLE with a harmonic potential, i.e. 
for any a < ur there always exist a specific fi/j (which 
depends on 7 and co) for which the system is in resonance. 
The second finding is that above aR we have a well de- 
fined boundary between the phase for which the reso- 
nance exist - a "Resonance Phase" , and a phase where 
are no resonance - "No-Resonance Phase" (see Fig. 111. 
For a < aR "No-Resonance-Phase" does not exists. The 
boundary is given by the following relation 



5(a), 



(81) 



and the phase diagram is presented in Fig. |ll[ The func- 
tion g{a) is found analytically (see Appendix IX I 



g{a) = c{ay 



/c(a)^ - c(a) cos (^)(l-ff) + f 



c(a) cos I 



TTOL ^ 
2 / 



(82) 
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where c(a) is given by Eq. ( 108 1. As is seen from Fig 
for u;/^^^ > g{a) a "Resonance Phase" is obtained and 
"No- Resonance Phase" for uj/'y^^ < g{a). In the limit 



12 



of a ^ 1 the boundary goes to the expected value for 
damped oscillator l/\/2- Near the critical point an the 
g{cx) drops to zero as a power-low with exponent 1/2 



aR. 



(83) 



The existence of the same critical for a free and a 
harmonically bounded particle is easily understood from 



the phase diagram in Fig. 11 Choosing the straight line 
w = 0, which represent the free particle, and going along 
this line when starting at the "Resonance Phase" we will 
cross to the "No- Resonance Phase" exactly for a = aR. 
Generally speaking the same critical or will be obtained 
for any kind of external force because it is determined by 
the internal properties of the surrounding medium rep- 
resented by the friction part in FLE. The phase diagram 



in Fig. ITT] has much in common with the phase diagram 
in Fig. [7j which is quite reasonable because of the strong 
connection between non-monotonic decay of the correla- 
tion Cx{t) and existence of a resonance. A presence of a 
resonance for small enough a emphasize the previously 
obtained result of non-existence of the over-damped limit 
for such a. Those properties are due to the same cage 
effect that we already discussed. 



C. Complex Susceptibility 



Eq. ( 78 1 for the complex susceptibility can be written 



in the following form 



{uj^ - n'^) + -fn°' cos (f) . 7f]«sin(^) (84) 

(^2 _ ^^2)2 + ^2fl2a ^ 27(^2 _ f^2)f^a ^OS (^) \uJ^ - fl^)^ + -/^fl^'^ + 2j{u;^ - n^)n°' COS (^) ' 



the real and the imaginary parts of the susceptibility 
are experimentally measured quantities for many systems 
and so it is interesting to explore their behavior for the 
FLE. In this subsection we will explore the behavior of 
the imaginary part x"{^) which is also called "the loss". 
From Fig. 12 we observe interesting behavior of x"(^) for 
different a's, not only one peak is present as is expected 
for the normal oscillator, but we observe a double peak 
phenomena for some a's. The double peak phenomena of 
"the loss" for super-cooled liquids and protein solutions is 
a well known phenomena [47", '48', ^49] and usually treated 
by the means of mode-coupling theory ,48] . We define 
two phases for the behavior of x"(ri) and in the follow- 
ing find the phase diagram for x"(ri). The first phase is 
the phase where x"(^) has only one peak - "One- Peak" 
phase and a "Double-Peak" phase, where a double peak 
phenomena is observed. 

In Appendix pC] we explored dx"{^)/dfl and searched 
for the boundaries between the "One-Peak" and "Double- 
Peak" phases. The result is presented in Fig. [TSj The 
boundaries between the phases are given by analytical 
functions gi{a) and 52(a) which are only a dependent. 
Two critical a are found for such a phase diagram, the 



first one 



0.527... for which the boundary 51(a) 



drops to zero, and a second one a^^ = 0.707... for which 



gi{a) and 52(a) coincide (see Fig. 12 1. We also must note 



that near a^-^ 5i(a) behaves like 51(a) oc (a 



1/2 



a behavior which was also observed for g{a) and Km (a) 
near the corresponding critical points. 

The double-peak phenomena in our case is explained 
in the same sense as was explained the existence of reso- 
nance for small enough a, and the disappearance of the 




FIG. 12: The imaginary part of the complex susceptibility 
for various a. A double peak phenomena is observed, (a) 
a = 0.66, 7 = 1.8 and lo = 0.7. (b) a = 0.5, 7 = 10 and 
uj = 2. (c) a = 0.63, 7 = 10 and = 2. (d) a = 0.8, 7 = 10 
and Lu = 2. 



monotonic decay phase for the correlation function. The 
reason is the same, the friction becomes more of an elas- 
tic force for such small a embedding oscillations in the 
system. Such a claim is emphasized by the Cole-Cole 
plots of the complex susceptibility, presented in Fig. [14] 
For a = 0.8 the behavior is very much as for a normal 
damped oscillator, a Debye mod el |16j . for small enough 
uj/^^^^^~"'> (large friction) - Fig.ll4l(c), i.e. a monotonic 
behavior for the relaxation, and a Van Vleck-Wcisskopf- 
Frohlich type [IW. for a large value of a;/7^/*^^~"' (small 
friction)- Fig. 14 (d), i.e. a oscillating behavior of the re- 
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04 0.6 



0.05 0.05 



FIG. 13: Phase diagram of the imaginary part of complex 
susceptibihty. Two phases are found, a phase with a presence 
of one maxima for x"(^) " "One-Peak" phase, or a phase 
with a presence of two maxima - "Double-Peak" phase. For 
a > = 0.707... only one phase present. 



laxation. These two prototypes of normal complex sus- 
ceptibility correspond to presence of a single character- 
istic frequency in the system, or a single time scale if we 
are concerned with correlations. For small a, like a = 0.1 



in Fig. 14 'a), we see a coexistence of these two types of 

'a) cor- 
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normal susceptibilities. The right side of Fig. 
responds to a Debye type - a monotonically decaying pro- 
cess and on the left side a Van Vleck-Weisskopf-Frohlich 
type wich shows highly oscillating behavior even in the 
case when 7 and uj are the same as for Fig. 14 (c). Effec- 
tively for small a we have two characteristic frequencies 
in the system, the lowest is responsible for the mono- 
tonic decay and the high frequency is an oscillating pro- 
cess, for intermediate a we have some mixed behavior 
- Fi g. [T4l (b). This oscillating behavior that is seen in 
Fig. |14[(a) is the manifestation of the cage effect which 
we already explained as the rattling motion of the sur- 
rounding particles and is presented in FLE due to the 
friction force. 



VII. SUMMARY 

The Fractional Langevin Equation (FLE) with power- 
law memory kernel and < a < 1 is a stochastic frame- 
work describing anomalous sub-diffusive behavior. This 
equation may be expressed in terms of fractional deriva- 
tive and so provides an example of a physical phenomena 
where non- integer calculus plays a central role. The so- 
lution of a fractional-differential equation describing the 
correlation function was presented in terms of roots of 
regular polynomials. It was shown that for a ^ 1 there 
is no unique way to define an over-damped or under- 
damped motion. Three definitions were proposed for the 
frequency of transition, i.e. ujc, and tOz- We ob- 
served an existence of a phase transition for a critical 
a = ac ^ 0.402, where for a < ac Cx{t) does not decay 



FIG. 14: The Cole-Cole plots of the complex susceptibility, 
as a function of {^) ~ 0.1, 7 = 10 and uj = \. 

(b) a = 0.5, 7 = 10 and (.J = 0.5. (c) a = 0.8, 7 = 10 and 
1^ = 1. (d) a = 0.8, 7 = 1 and a; = 10. Panel (c) and (d) 
show two typical normal behaviors, for small a the panels (a) 
and (b) show a behavior which is a mixture of the two normal 
typical behaviors. 



monotonically for any > 0. Physically it is explained 
using a cage effect : a rattling motion of a particle in the 
cage formed by the surrounding particles. A response to 
a time dependent field in terms of complex susceptibil- 
ity x(^) also was calculated and similar critical a's were 
found. Particularly for a < an = 0.441... the system 
will always be in a resonance with the external field for 
particular f2 and any 7 and lo, even in the case of free 
particle {uj = 0). For the "loss" - the imaginary part of 
the complex susceptibility, x"(ri), two phases were de- 
fined: (i) "One-Peak" phase were the complex suscepti- 
bility obtains only one maximum as in a regular case, 
(ii) "Double-Peak" phase, were the complex susceptibil- 
ity obtains two maxima, a phase diagram was presented. 



and = 0.707. 



Two critical exponents a^-^ = 0.527. 
were found for x"(f2), exponents which define the bound- 
aries of the phase diagram. In conclusion, critical expo- 
nents like ac, ctfi, a^^ and a^^, mark sharp transitions in 
the behaviors of systems with fractional dynamics. Thus, 
these critical exponents are clearly important and general 
in the description of anomalous kinetics. 
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VIII. APPENDIX A: THE SOLUTION FOR 
NON-DISTINCT ZEROS OF P(s) 



In this Appendix we derive the solution of Eq. ( 12 1 for 
the case when two zeros of P{s) coincide. This means 
that P(s) has 2q — 2 distinct zeros of order 1 and one 
zero of order 2. Namely, at the critical point ujc only two 
Qk coincide and we present here a method of solution for 

LU — UJc- 

Starting with Eq. (12 1, we write the partial fraction 
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expansion in the following way 

2g-l 



1 



E 



A, 



A 



P{s) ^ (s- 029-1)' 



(85) 



where Ofc are the zeros of P{s) and we assign 02^-1 to be 
the zero of the second order. and A are given by 



Using the relation 



2g-l 



^ al'Ak + mal'-^A = ?7i = 0, 1, . . . , 2q~2 (89) 



k=l 



Ak = 



dP{s) 
ds 



l<fc<2g-l, (86) 



one finds that 



and 



A2q-1 



A = 



2(j-2 
k=l 



1 



J |s=a2„-i 

as s — a2q-i — 



(87) 



2g-l 

E 



k=l 



S - flfe 



P(,s) 

Hence using Eqs. ( |12|18|90 l 



S - a2q-l [S - a2q-iy 

TO = 0, 1, . . . ,2g - 1. 



(90) 



2g-l(j-l /2g-l 

c'^(^)=EE E 

m=0 i=0 \ fc=l 



_i iT^o,2„_iA]3„ij _j_ a2„_iAB,yij 

S 1 + S t + - — r^S <J 



S - a2q-l 



and it is only left to perform an Inverse Laplace Transform of — , 

S9 (s - a2g_i)2 



{s-a2q-lY J 

using convolution theorem 



S' (S - a2q-lY 



or using Mittag-Leffler function 



1 



Sl{s~ a2q-lY 



7 7(-,a2q_it) rrr7(- + l,a2q-it) 

r(|)a|,_i « r(i)a|;_i 1 



t'i^^E. J {a2q-lt) - -t^^^E. (a2g-lt) ■ 



Finally using Eqs. ( |21|91|92 l 



m=0 j=0 



■2g-l 



^ a" Mfce'^'=*7(^,afet) + a",_fie'^^'-^* ( (< + ma^,-i)7(-, «2g-it) - <-i7( 
fc=i \ 9 9 



J 



(91) 



(92) 



(93) 



or using Eq. ( 93 1 



l,a2q-ii) 
(94) 



29-1 q-l 



m— — 



■29-1 



E <^kE, ,^^{akt) + a%_^A {t + ma2^_^)E^ (02,-1*) - -tE^ 2+^- («2, 



q-l 



t) 



r 



(95) 



A final remark: one can show that for our case of inte- IX. APPENDIX B: EXPLORATION OF EQ. ( [80| ) 
ger g > p > 0, third and higher order zeros of P(s) don't 
exist. 



In this Appendix we prove the existence of for FLE 
with a harmonic force and derive the equation for g{a) 
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given by Eq. Solving Eq. ^ one gets for 7/f^|"" 

' ^([2-«.]cos(f); 



2-a 



^2a V " "yPcos (™))' + Say - 'Zi±(y), 



(96) 



and y — ( — 1 ) > — 1. Writing the left hand side of 



Eq. ( 96 ) in terms of y 



7 



7 



ni 



q2{y)- 



(97) 



We see that for the extrema points of i?(il) the functions 
qi±{y) and (72(2/) cross each other (see Fig. 15 1. While 



92(2/) is a monotonic increasing function starting from 
zero for y = —1 and growing as y^~^ for large y, qi±{y) 
constructs two branches where qi+{y) is the upper branch 
and qi^{y) is the lower branch and for some point y^ 



qi-iVa) = qi+iVa) 



(2-aj/„)cos(^) 



2a 



(98) 



If j/q < —1 then 92(2/) crosses qi±{y) no mater what 
the parameters 7,0; and a are, because in that case for 
y = —\ qi±{y) > and 92(2/) = and a resonance is 
always obtained. The point ya is derived from Eq. (96) 
and determined by the following relation 



(2-ay„)2cos2 ( 



7ra\ 



+ 8aya = 0. 



Solving Eq. ( 99 1 in terms of ya one finds 



a cos- 



/ 7ra \ 
\ 2 J 



TTa\ 



(99) 



(100) 



where we took the — sign because y > — 1. For < a < 1 



Eq. ( 100 1 is an increasing function of a and so we have 
a critical a, aj^ for y^ = —1. Eq. (99) with y^. 



-1 is 



we have shown the existence of an for the harmonically 
bounded particle. 

We now argue that the boundary between "Resonance- 
Phase" and "No-Resonance Phase" is given by the follow- 
ing relation 



where g(a) is only a dependent and equals zero for 
d < Q-R- One readily sees from Eq. ( [96| that the 
large y behavior of qi±{y) is proportional to ±y^, 
where (72(2/) behaves like y^^ ' 
also we note that '^'^^^''^^ 



dy 



and ^ 



for large y (Eq. (|97|), 
are monotonically 

decaying functions of y. As a result we have four options 
for the scenario of q2 (y) crossing qi± (y) (see also Fig. 15 1 : 



(i) 92(2/1) = '7i-(2/l) and 92(2/2) = 9i+(2/2) for yl < y2, 



-1 -0.5 D.5 1 

m 
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V 
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y 



FIG. 15: The four scenarios for qi±{y) and q2{y) crossings. 
Existence of a crossings correspond to "Resonance-Phase" 
and no crossings correspond to "No- Resonance Phase" . (a) 
a = 0.441, 7 = 2 and lu = 1, for this case there will always be 2 
crossings, a < an. (b) a — 0.6, 7 = 1.75 and uj — 1, this panel 
corresponds to a "Resonance-Phase" (c) a = 0.5, 7 = 4.4 and 
ijj = \, this panel describes a situation on the boundary be- 
tween "Resonance-Phase" and "No- Resonance Phase", (d) 
a = 0.5, 7 = 10 and tj = 1, a "No-Resonance Phase". 



(ii) '72(2/1) = 9i+(2/l) and 92(2/2) = '?i+(2/2) for yl < y2, 

(iii) 92(2/) = q\+{.y) for a single y 

(iv) (72(2/) 7^ 91+ (2/) and 172(2/) 7^ 9i-(2/) for any y. 

When there are two crossings then the one with 
the larger y corresponds to minimum and the smaller y 
corresponds to maximum and belongs to the "Resonance- 
Phase" . When there are no crossings then we are in the 
"No- Resonance Phase". The scenario (iii) corresponds 
exactly to the boundary between the two phases. In 
order to find the boundary, two conditions must be 
fulfilled 



92(2/1) = 9i+(2/l) 



exactly Eq. ( 77 ) which defines the equation for and so ^^(j 



dq2{y) 
dy 



la=yi 



dqi+{y) 

dy 



(102) 



(103) 



as illustrated in panel (c) of Fig. 15 Starting form 
Eq. (102) we compare the left-hand side to some con- 



(101) 

stant c and using Eq. (|97|) we find 



- 1. 



(104) 



Comparing the right hand-side of Eq. (102) to the same 



c and using Eq. (104) we arrive to the following relation 



ccos 



(¥)(l+^) + ^ 



;cos {f) 



(105) 
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Next performing the derivation in Eq. ( 103 1 and using tion for c 



Eq. ( |104| we find 

/ \ 2 



(2-a)c(2c-(l + f cos(^)) _ . 



.7 = 



4 



4ccos(^) 



ac cos I 



cos (^) + (2a - - a{2 + a) cos (Tra)) 0^+ 
12acos(^)c-8 = 0. 

(107) 

(106) Eq. (1071 has three different solutions where only one is 



comparison of Eq. ( 106 1 and Eq. ( 105 1 supplies an equa- real for < a < 1, we will call it c(a) and 



c{a) 
1 

3^ 



sec (-y-) {—2a + ^c? + 2q;cos {-ko) + o? cos (vra) 



{2^c?^m^ (-4 + 20a-7a2 + (2 + a)2cos(7ra))' 

[-80a^ + 312a^- 384a^ + 1520;*^ - io? (-40 + 132a - 126a^ + 43a^) cos(7ra) - 48a^ cos(27ra) + 72a'' cos(27ra) 
— 24a^ cos(27ra) + 8a'^ cos(37ra) + 12a^ cos(37ra) + 6a^ cos(37ra) + a® cos(37ra) 



24\/6j-a9 cos2 



("^) (^^ ~ '^^'^^ ^ '^^^^'^ ^ ^^"^ + (2 + a)2(-8 + 5a) cos(7ra) sin* 



([-80a^ + 312a'* - 384a^ + 152a'5 - io? (-40 + 132a - 126a2 + 43a^) cos(7ra) - 48a^ cos(27ra) + 72a"* cos(27ra) 
-24a^ cos(27ra) + 8a'^ cos(37ra) + 12a'* cos(37ra) + 6a^ cos(37ra) + a^ cos(37ra) 



24\/6w'-a9 cos2 (^32 - 204a + 204a2 - 59a3 + (2 + a)2(-8 + 5a) cos(7ra) sin" 



7ra\ 
~2) 



1 \ n 

3 



(108) 



We thus justified the use of Eq. (101 1 and g(a) is given 
by Eq. (l82|. 



X. APPENDIX C: EXPLORATION OF dx\^)l<Kl 



We start with the exploration of d'yl' (^) j d£L, where 
X"(r!) is given by Eq. (Istl), 



dfl 



a+3 



sm 



\ 2 ) 



f]2 



7ra\ 



- -4 



/TTa 



(109) 



and we easily see that in order to '^^J^^^ = 0, the follow- 
ing condition must be fulfilled 

(110) 

where y = ^fj2— 1^ > —1- The left hand side of 



Eq. (110 1 is a second order polynomial in terms of 



o2- 



which is easily solved 

( ' 



.^2 



4cos(^' 



2a 



± — a/16cos2 (Jy) +4a(4y + ay2). 



1 

'2^ 



(111) 

The right hand side of Eq. (Ill I we will call qi±{y) and 
the left hand side 52(2/), 



92 (y) 



7 



,2-a 



(y+1) 



1-% 



(112) 
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The crossings of qi±{y) and q2{y) determine the extrema 
points of x"(ri) and using the fact that the for y oo 
Q2{y) y^^T and qi±{y) oa y we have six different 
scenarios for the crossing of 92(2/) and qi±{y). 

(i) 92(2/) = qi-{y) for a single y and there are 
no other crossings, 

(ii) 92 (yi) = 91- (2/1) and ^2(^2) = 91+ (2/2) for 
2/1 < 2/2, 

(iii) 92(2/1) = <7i-(2/i), 92(2/2) = 9i+(2/2) and 
92(2/3) = 9i+(2/3) for 2/1 < 2/2 < 2/3, 

(iv) 92(2/1) = 9i+(2/i), 92(2/2) = 9i+(2/2) and 
92(2/3) = 9i+(2/3) for 2/1 < 2/2 < 2/3, 

(v) 92(2/1) = 9i+(2/i) and 52(2/2) = 9i+(2/2) for 
2/1 < 2/2, 

(vi) 92(2/) = 9i+(2/) for a single y and there are no 
other crossings. 

Generally if there is only one crossing, scenario (i) 
and (vi), the meaning is that x"(ri) will have only one 
maximum and on the contrary when there are three 
crossings, scenario (iii) and (iv), there are two maxi- 
mums and one minimum for 'x"{^). These correspond 
to two different phases the "One-Peak" phase and the 
"Double-Peak" phase where the scenarios (ii) and (v) are 
the boundaries between these phases. We are interested 
in finding these boundaries, where for scenario (ii) and 
(v) two conditions must be fulfilled 



92(2/1) = 9i+(2/i) 



(113) 



Next performing the derivation in Eq. (1141 and using 
Eq. (fUSf we find 



+ + 2 (^) " 4COS (^))) c-^ 

comparison of Eq. (116) and Eq. (1171 supplies an equa- 
tion for c 



2=4 



16 COS (^) {2 + a)c^ 
16 (-4 + 8a-a^ + cos^ (^) (2 + a)^) 5^ (Ug) 



cos(^f)(6-a)c+||(4-a) = 0. 



Eq. (118 



is a forth order polynomial and could pe solved 
by standard methods or using Mathematica. It has 4 
different solutions while two of the solutions have non- 
zero Imaginary parts for any < a < 1, while the other 
two solutions have no Imaginary part for a < a^^ ~ 
0.70776. Lets call these solutions Ci(a) and £2(0), the 
non-zero Imaginary part for a > of both c\ (a) and 
£2(0;) means that only scenario (vi) is applicable for such 
a's and we are in the "One-Peak" phase. The boundaries 
between the two phases are given by 



and 



dhjy) _ dqi+{y) 

dy \y = y. 



dy 



\y=yi 



(114) 



Starting from Eq. ( 113 1 we compare the left-hand side to 



some constant c and using Eq. (1121 we find 

, ,2 



2/ 



1. 



(115) 



1' 



\ 



5i,2(a) 




"^'"'l ' '' 51.2(a) - 5?,2(a) )ci,|" 



(119) 

where the subscript 1 is for gi{a) the lower bound in 
Fig. 13 and subscript 2 is the upper bound 32(0) in 
Fig. 13 For gi{a) there is also another interesting point 



Comparing the right hand side of Eq. (1131 to the same a^i = 0.527031 which satisfies the following relation 



c and using Eq. (115) we arrive to the following relation 




-4 



4 cos(^ ^) ^ 



a^^-4a,i+4cos2('^)=0. 



(120) 



(116) for a < 51(0) — 0. 
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